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We consider a time-dependent non-linear Schrodinger equation in one dimension (ID) with a 
fifth-order interaction term and external harmonic confinement, as a model for both (i) a Bose 
gas with hard-core contact interactions in local-density approximation, and (ii) a spin-polarized 
Fermi gas in the coUisional regime. We evaluate analytically in the Thomas-Fermi limit the density 
fluctuation profiles and the collective excitation frequencies, and compare the results for the low- 
lying modes with those obtained from the numerical solution of the Schrodinger equation. We find 
that the excitation frequencies are multiples of the harmonic-trap frequency even in the strong- 
coupling Thomas-Fermi regime. This result shows that the hydrodynamic and the coUisionless 
collective spectra coincide in the harmonically confined ID Fermi gas, as they do for sound waves 
in its homogeneous analogue. It also shows that in this case the local-density theory reproduces the 
exact collective spectrum of the hard-core Bose gas under harmonic confinement. 
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I. INTRODUCTION 



The study of atomic gases in condensed quantum states inside magnetic or optical traps has received an enormous 
impulse from the achievement of Bose-Einstein condensation in vapors of bosonic alkali atoms and of atomic hydrogen 
. Similar techniques of trapping and cooling are being used to bring gases of fermionic alkali atoms into the quantum 
degeneracy regime |2j. The observation of small-amplitude shape-deformation modes in Bose-Einstein condensates 
inside three-dimensional (3D) anisotropic traps at extremely low temperatures js] has proved to be an important 
method of diagnostics and has provided a crucial test of the mean-field theory based on the time-dependent Gross- 
Pitaevskii equation |Q. No such experiments have as yet been observed in Fermi gases, but predictions are available 
on their dynamical spectra in both the coUisional and the coUisionless regime For both bosons and fcrmions in 
3D the strong-coupling Thomas-Fermi limit is characterized by a different spectrum from the weak-coupling limit. 

Attention has also been moving towards atomic gases in restricted geometries, where new frontiers are met in the 
realization of atom lasers [^|J^] and of thin atom waveguides The rich phase diagram of a quasi-lD Bose gas may 
become accessible to observation through tuning of the interactions |^ . Essential modifications to the Gross-Pitaevskii 
equation are needed in D< 2, since the T-matrix vanishes already in D=2. The dilute gas limit can still imply strong 
coupling in ID, since collisions are unavoidable in this reduced dimensionality. 

At low temperature and density the dynamics of a Bose gas with repulsive interactions inside a very thin waveguide 
approaches that of a ID fluid of impenetrable point-like bosons As demonstrated by Girardeau this so- 
called Tonks gas pT| has the same spatial profiles as those of a spin-polarized Fermi gas, since the exact many-boson 
wavefunction is related to that of non-interacting spinless fermions by b{xit..,xn) = |^'F(a^i, ••■7 a;jv)|. More 
generally, the Fermi-Bose mapping theorem implies that the dynamic structure factor of the Tonks gas is the same as 
that of the corresponding ideal Fermi gas, although their momentum distributions are very different. In particular, 
in the same work Girardeau [|lO| also showed that long-wavelength sound waves propagate in the homogeneous Tonks 
gas with velocity v = nhn/m equal to the Fermi velocity vp in the homogeneous ideal Fermi gas at the same density 
n and atom mass m. This is, in fact, the speed of propagation for both zero (coUisionless) sound and hydrodynamic 
sound in the ideal ID Fermi gas. 

The dynamic structure factor of the spin-polarized ID Fermi gas under harmonic confinement has been evaluated 
exactly in the coUisionless regime and related to that of the homogeneous gas by means of a local-density formula 
| p2| . The collective excitation frequencies in this regime are integer multiples of the harmonic trap frequency and this 
spectrum is, therefore, also that of the corresponding Tonks gas. In the present work we evaluate the dynamics of 
the same Fermi system in the coUisional regime described by linearized hydrodynamic equations. We shall show how 
these equations for the Fermi gas can be derived from the non-linear Schrodinger equation proposed by Kolomeisky 
et al. [l3| for the dynamics of the mesoscopic wave function of the corresponding Tonks gas and how the equality of 
the spectra of the two fluids remains valid within the dynamical regime of present interest. 

After a presentation of the model in Section II, we evaluate its dynamical density fluctuations and collective 
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frequency spectrum in the following Sections. We give in Section III an analytic solution of the equations of motion 
for the density fluctuations, using a Thomas-Fermi approximation (TFA) extended, however, to consider solutions 
external to the classical radius. In Section IV we treat the same problem by means of a numerical simulation mimicking 
an experimentally realizable method for measuring the frequencies of the gas. That is, we excite a collective mode of 
the atomic cloud by applying an external time-dependent potential of chosen symmetry and probe it by observing the 
evolution in time of the density profile once the perturbation is turned off. Finally, Section V gives a brief summary 
of our main results and some concluding remarks. 



II. THE MODEL 



According to Kolomeisky et al. ||T^, a long- wavelength approach to the hard-core ID Bose gas with repulsive 
point-like interactions in the dilute regime can be based on the non-linear Schrodinger equation 



ihdt^(x,t) 



Zm Zm 



<f(x,t). (1) 



Here, V{x,t) — Vextix) + Up{x,t) is a time-dependent external potential, which includes the harmonic confinement 
Vext{x) = mijj\^x^ jl and a periodic perturbation Up{x,t). The wave function $(a;,i) is normalized to the number 
N of particles in the trap. Equation takes account of the correct density dependence of the ground-state energy 
density for the ID Bose gas in the strong-coupling limit within a local-density approximation. 

It was shown by Girardeau and Wright that Eq. (|l|) describes well the dynamics of the Bose gas, although it 
overestimates its coherence in interference patterns at small numbers of particles. Equation (|l|) can be cast in the 
form of Landau's hydrodynamic equations for a superfluid by setting ^{x,t) — \Jn{x^t) exp[i(^(a;, i)]. This yields 

dtn{x, t) = -dx [n{x, t)v{x, t)] (2) 

and 
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mdtv{x, t) 



liioc{x,t) + ^mv^{x,t) 



(3) 



where the velocity field is defined as u(x,t) — {h/m)dx4>{x,t) and the local chemical potential is given by 



K TT h 

tiioc{x,t) = ~- ^==dl^nix,t) + Vix,t) + -—n^ix,t) . (4) 

2m^n{x,t) zm 

The TFA is obtained by neglecting the first term in the RHS of Eq. (|^) . Notice that the ratio of the kinetic energy 
to the interaction energy is proportional to iV"^. The fact that in the TFA the local chemical potential scales with 
density as 'k^W't? j^ra^ as for a non-interacting spin-polarized Fermi gas in the local-density approximation, reflects 
the property that the hard-core Bose gas has the same spatial profiles as those of the Fermi gas. 

The boson-fcrmion mapping has also been extended to time-dependent phenomena pO|Jl^ . Therefore, using the 
mapping in the inverse direction we expect that Eq. describes approximately a non-interacting Fermi gas within 
a local-density approach. Indeed, it is easily shown that upon neglecting the kinetic energy term in fiiodx^t) Eqs. (j^) 
and (|^) yield the equation of motion for the density fluctuations of the Fermi gas in the hydrodynamic regime. The 
latter is given by 

mdtn{x, t) = dllV{x, t) + dx [n{x, t)dxV{x, t)] (5) 
where the momentum flux density Il{x,t) depends locally on the particle density through the relation 

n(x,t) = h'^TT^n^{x,t)/STn + mn{x,t)v^{x,t)/2 . (6) 

Equations (||) and (|^) are obtained quite straightforwardly by combining Eqs. (|^)-(^) in the TFA. 

Having established that Eq. (|l|) describes in the TFA limit the dynamics of both hard-core bosons in the strong- 
coupling regime and non-interacting spinless fermions in the hydrodynamic regime, we turn to its equilibrium solution, 
which is a necessary preliminary step for obtaining the collective modes. 

The ground-state density no(a;) is obtained from Eq. (^) by setting ^ioc{x) = with v = Q and Up — 0. This yields 
the following differential equation for the equilibrium density: 

+ -^'^o^^(^) + ["^extix) - n]^/na{x) = . (7) 



2 



It is important to remark that Eq. (|^) does not yield the exact equilibrium density profile of the gas even if one 
includes the kinetic energy term {h"^ /2m)d^ ^ uq (x) . The exact profile satisfies instead the third-order differential 
equation |16 



-—dlno{x) - -no{x)dxVext{x) + [Vext{x) - fj]dxno{x) = . (8) 
8m 2 

In the case of = 100 particles an accurate comparison between the solutions of Eq. (Q) (obtained by the numerical 
method described in Sec. |^ below) and of Eq. (g|) (using the results of Ref. |l^) is given in Fig. |l|. It is seen that 
although the overall shape of the two curves is quite similar, Eq. (^ misses the shell structure of the exact profile and 
has an incorrect behavior in the tails (sec the insets in Fig. |^) (see also Ref. |p^). 

The above comparison illustrates the limits of validity of Eq. (|^) . We shall sec in the following that the spill-out 
of the particle density beyond the classical radius plays a crucial role in determining the dynamics of the atomic 
cloud. We may also remark that, whereas Eq. (|^) is constructed from an approximation for the kinetic energy density 
functional, the exact functional is in fact explicitly known for this system [l^ ]. However, it involves a complicated 
non-local function of the particle density. 



III. COLLECTIVE MODES IN THE THOMAS-FERMI LIMIT 



The linearized equation of motion for small-amplitude density fluctuations is 

mdf5n{x,t) = dlSUix,t) + dx[Snix,t)dxVext{x)] , (9) 

where 5n(x,t) = n{x,t) — uq^x) and we have assumed resonance conditions, i.e. Up{x,t) = 0. In the homogeneous 
gas Vext{x) — and Eq. (|^) yields an acoustic dispersion relation with the propagation velocity of hydrodynamic 
sound. It is a peculiarity of the ideal ID Fermi gas that this coincides with the velocity of coUisionless sound. 
In the strict TEA limit the equilibrium density profile has finite extension, as it is given by 



H ' (x) = ^0{xl - x')[^l - v,xt{x)\^'^ (10) 

where the chemical potential ix = HuJhoN is obtained from normalization of Eq. ( [lO| ) to the number N of particles 
in the trap. We define the classical turning point as the Fermi radius Xp — \/2Naho with a^o — ^/hjmujho. 
Correspondingly the TEA expression for the fluctuations in the momentum flux density is 6IlTFA{x,t) = 29 {Xp — 
x^)[/i — Vext{x)]Sn{x,t). As will be apparent from the arguments given below, the spectrum of collective excitations 
is obtained as a continuum if this expression is used in Eq. ^ . 

A discrete spectrum is instead obtained if the spill-out outside the classical radius is taken into account. To do this, 
we adopt the approximate relation 

dll{x, t) ~ 2^6n{x, t) = 2[^l - Vext{x)]5n{x, t) (11) 
on 

where t{rL) is the kinetic energy density and the second equality follows from using the Euler equation for Density 
Functional Theory. 

By substituting Eq. ( pT| ) into Eq. (^ and performing a Fourier transform with respect to time, we rewrite Eq. (|^) 
as the eigenvalue equation 

(1 - z^)dl5n{z, oj) - ZzdM{z, oj) + [{uj/ujhof - l]5n{z, w) = , (12) 

where z = xjXp is a rescaled position variable. It can be checked by direct substitution that for \z\ < 1 two 
independent solutions of this second-order differential equation are 

(1) cos[(cj/wfto)arccosz] 
SKn (2. ^) = yj— ^ (13) 



and 



Sn^^iz,co)^'^^^^M!^^^^S^ (14) 
Vl — 
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for any value of the ratio uj /iOho- In the domain \z\ > 1 we instead have 

dnout(z,uj}^ , (15) 

yz — 1 

having imposed that the density fluctuations vanish for \z\ oo. Finally, the dispersion relation and the shape of the 
density fluctuation modes inside the classical radius are obtained by imposing continuity of -^/jl — z'^\Sn{z, lu) across 
the Fermi radius. This condition selects the integer values of the frequency ratio, 

t^/w/io = n (16) 
and the form ( |l3| ) of the inner solution. This can also be written as 

X (1)/ \ T.n{z) 

6nl^'{z,nuJho) ^ ^^^== (17) 

where the functions T„(z) are the Chebyshev orthogonal polynomials of the flrst kind iQ. The (integrable) divergence 
in Eq. ( |l7|) at the Fermi radius is a consequence of the linearized TFA in low dimensionality and is evidently 
inconsistent with the linearization of the hydrodynamic equations. We shall return on this point in the next Section. 

In the case n = 1 we correctly recover from Eq. ( p^ the frequency of the sloshing mode, in agreement with the 
generalized Kohn theorem ||l^ . More generally, we find that the TFA spectrum coincides with the spectrum of the 
Fermi gas in the collisionless regime, as evaluated in Ref. jl^. Although this result for the trapped gas has been 
reached by an approximate and to some extent ad hoc argument, it will be verifled immediately below by direct 
numerical solution of Eq. (m) . 



IV. NUMERICAL SIMULATION 



The divergence of the density fluctuations at the classical boundary of the cloud that result from the extended TFA 
are incompatible with a linearization of the equations of motion. We have therefore carried out a numerical solution of 
Eq. (H), which reproduces an experimental procedure employed to excite collective modes [1). Namely, we first apply 
a time-dependent perturbing field of given symmetry for variable amounts of time corresponding to several periods 
of the expected excitation, and then monitor the evolution of the cloud by recording a series of "pictures" once the 
perturbation is turned off. From the evolution of the density profiles we obtain the frequency and the shape of the 
density fluctuation modes. 

The numerical simulations use an explicit time-marching algorithm pO| | for propagating the state of the cloud 
in both imaginary and real time, yielding the ground-state proflle no (^j and the subsequent dynamical behavior 
n{x,t). Due to the size of the nonlinear term in Eq. (|l]), special attention has been given to the stability of the 
algorithm. The perturbing potential has the form Up{x,t) — Up{x) cos{uJit) and the drive frequency uji is tuned 
around the expected resonances. As to the spatial variation of the perturbation, we have adopted the TFA form 
Up(x) cx a/I — [x/XpYSn]^ (x, oji) in order to excite well defined modes. This is approximately equivalent to imposing 
orthogonality between the drive field and the other density fluctuations, in view of the orthogonality of the Chebyshev 
polynomials in the domain \z\ < 1. 

Figure ^ shows the spectra of quadrupole and hexapole collective modes for N = 100 particles as obtained from the 
Fourier transform of the time evolution of the density profile taken at a given point in the harmonic trap. The results 
agree with the TFA prediction that the main frequencies of these modes are w = 2 ujho and lj — i uto, respectively. 
Nonlinearity is apparent through the presence of several other resonances in each spectrum. We have checked that 
the spectra of the i*^ moments (x*) = (1/A^) / dx x^n{x, t) are peaked at the fundamental excitation frequency of each 
mode. In further tests we have (i) examined the behavior of the quadrupole mode frequency with varying N, finding 
only very small deviations within the error bars from the value u = 2ojho at low N\ and {ii) used non-integer values 
for the driving frequency, finding again only spectral lines at integer frequencies. 

Finally, in a regime of weak drive amplitude we have compared the density fluctuation profiles in Eqs. ( p^ and ( p^ ) 
with the output of the simulation. The comparison is illustrated in Fig. ^. The TFA density fluctuations, in addition 
to suffering from divergiencies at the classical radius, do not integrate to zero for the modes with n =even integer. 
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V. SUMMARY AND CONCLUSIONS 



In summary, we have shown by analytical and numerical arguments that the hydrodynamic frequency spectrum of 
the ID ideal Fermi gas under harmonic confinement is given by integer multiples of the trap frequency. This result 
extends to the trapped Fermi gas a property of the homogeneous ideal Fermi gas, in which the speed of sound is the 
same in the hydrodynamic regime as in the coUisionless regime. A more accurate treatment of the collisional regime 
beyond the TFA may usefully be developed in the future. 

A further conclusion can be drawn from the present solution of the non-linear Schrodinger equation describing a 
ID hard-core Bose gas. By virtue of the fermion-boson mapping theorem, the frequency spectrum associated with 
this local-density theory reproduces the exactly known spectrum of the Bose gas under harmonic confinement. 
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FIG. 1. Equilibrium density profile for A'^ = 100 fermions, in units of the harmonic oscillator length auo = \/ h/mtoho- The 
exact profile (solid line) is compared with that obtained from the non-linear Schrodinger equation (fi]) (dashed line) . The insets 
show enlargements of the region around j: = and of the tail region, in the same units. 
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FIG. 2. Spectra of density fluctuations obtained from the numerical simulation by applying a perturbing potential of drive 
frequency u — 2 ujho (left panel) and uj = ?> uiho (right panel) and of the appropriate symmetry as described in the text. The 
insets show the corresponding time evolution of the density profile n{x,t) (in units of aj^^, with t in units of uj^^) taken at a 
given point x in the trap. The spectra are obtained from the Fourier transform of n{x,t), extended over several periods for a 
better definition of the mode frequencies. 
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FIG. 3. Density fluctuation proflles 5n{x) for the low-lying modes (indicated by the mode numbers n = 1, 2 and 3 from 

top to bottom) at A'^ = 100 particles, in units of a'^^ as functions of position x in units of the Thomas- Fermi radius Xf- The 
results from the numerical simulation (solid lines), where 5n{x) = n{x,t) — no{x) and t is a given time, are compared with the 
analytic solutions in the Thomas- Fermi limit (dashed lines). 
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